Classical momentum diffusion in double-5-kicked particles 
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We investigate the classical chaotic diffusion of atoms subjected to pairs of closely spaced pulses 
('kicks') from standing waves of light (the 2(5-KP). Recent experimental studies with cold atoms 
implied an underlying classical diffusion of a type very different from the well-known paradigm of 
Hamiltonian chaos, the Standard Map. The kicks in each pair are separated by a small time interval 
e ^ 1, which together with the kick strength K, characterizes the transport. Phase space for the 25- 
KP is partitioned into momentum 'cells' partially separated by momentum-trapping regions where 
diffusion is slow. We present here an analytical derivation of the classical diffusion for a 2(5-KP 
including all important correlations which were used to analyze the experimental data. We find a 
new asymptotic {t oo) regime of 'hindered' diffusion: while for the Standard Map the diffusion 
rate, for K ^ 1, D ^ K'^ /2[1 — 2J2{K)...] oscillates about the uncorrelated rate Do = K'^ /2, we find 
analytically, that the 25-KP can equal, but never diffuses faster than, a random walk rate. We argue 
this is due to the destruction of the important classical 'accelerator modes' of the Standard Map. 
We analyze the experimental regime 0.1 < Ke < 1, where quantum localisation lengths L ~ fi""'^* 
are affected by fractal cell boundaries. We find an approximate asymptotic diffusion rate D oc K'^e, 
in correspondence to a D oc K"^ regime in the Standard Map associated with 'golden-ratio' cantori. 

PACS numbers: 05.60.-k, 05.45.Mt, 05.45.Ac, 32.80.Pj 
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I. INTRODUCTION 

The 'i5-kicked particle' (5-KP) is one of the best known 
examples of classical Hamiltonian chaos. A particle, 
or in an experiment usually a large ensemble of ultra- 
cold atoms, is periodically 'kicked' by a series of very 
short laser pulses forming standing waves of light. The 
effective potential takes a sinusoidal form V{x,t) — 
—K cosxJ2n S{i ~ NT). Here T is the kicking period, 
while K is the kick strength, related to the laser inten- 
sity. The classical dynamics for the d-KP are given by the 
textbook example of chaos which is the 'Standard Map' 
0. For large K, the dynamics are characterized by dif- 
fusion in momentum. To lowest order, this represents a 
random walk in momentum , hence < >~ Dgt, where 
Do ^ K'^/2. The quantum counterpart of the S-KP is the 
quantum kicked particle (QKP). It has also been exten- 
sively investigated in numerous theoretical (see eg_Q,|3) 
and experimental studies in Austin |^ , Auckland 0, Ox- 
ford 01, Lille 0, and Otago 0. 

However, a recent experimental and theoretical study 
[Toj of cold cesium atoms exposed to closely spaced pairs 
of pulses (the 2(5-kicked particle) showed chaotic classical 
diffusion quite different from all other previously studied 
(5-kicked systems. The two kicks in each pair are sepa- 
rated by a short time interval e << T. 

The cold cesium atoms are of course a realization of 
the quantum counterpart of the 2S-KP; like the single- 
kick QKP the experiment exhibits the quantum chaos 
phenomenon of dynamical localization 2, 3, 4] whereby 
the quantum diffusion is arrested at a characteristic 
timescale, the 'break-time', t* ~ K^/h^. The momen- 
tum distribution of the atomic cloud is 'frozen' for times 
t > t* with a momentum localization length < >~ L^. 
By adjusting e and hence the timescales of the diffusion 



correlations, relative to the break-time, it was found in 
|lOl | that the experiment probed distinct diffusive regimes 
characterized by different 'families' of long-ranged corre- 
lations. In the Standard Map, the main corrections to 
the uncorrelated random walk come from 2 and 3-kick 
correlations. For the 2d-KP there were also families of 
terms correlating all kicks. Hence these additional 26- 
KP corrections were termed 'global correlations'. 

A particularly interesting feature of the 2d-KP is the 
'cellular' structure of the classical phase space. This 
structure is analyzed in detail below, but the basic idea is 
illustrated in Fig^ An ensemble of particles, all initially 
with momentum p = po diffuse chaotically, but it is seen 
that the diffusion is hindered at 'trapping regions' ie re- 
gions with momentap ~ ±(2m+l)Tr/e and to = 0, 1, 2, ... 
For small Ke trajectories spend most of their time stuck 
in the trapping regions. Once the particles have escaped, 
the timescale taken to diffuse freely over the remainder 
of the cell is negligible. 

The experiments in [lol | showed that for short and in- 
termediate times (~ 100 kicks for the experimental pa- 
rameters) the diffusion rates depended strongly on po and 
t. Here we examine also for the first time the asymp- 
totic regime (t — > oo) where the momentum spread of 
the atomic cloud is large compared with a single cell. At 
very long times momentum-dependent correlations de- 
cay to zero, leaving linear diffusion corrections leading to 
< p^ >= Doot- Fig [3 provides a summary of the short- 
time and the asymptotic momentum diffusion regimes in 
the 25-KP. 

The regime 0.1 < Ke < 1 is of special interest because 
it corresponds roughly to the experimental parameters 
and because a new study of the quantum 2S-KP [Til 
found a scaling behavior quite different from the usual 
QKP for this regime: while for the QKP the localization 
length L - /i-i, for the 25-KP, L - h'^ '^^. To date this 
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FIG. 1: (top) Surface of section plot for the 25-KP, K = 7, 
e = 0.05, T = T — e = 1.95, for atoms all with initial mo- 
mentum po = 0. The cellular structure is evident: momen- 
tum space is divided into regions of fast momentum diffusion 
separated by porous boundaries, ie narrow trapping regions 
where classical trajectories 'stick' for relatively long periods. 
The trapping regions are at momenta p ~ ±(2m-|-l)7r/e where 
m = 0, 1, 2, ...(bottom) A typical trajectory of the 25-KP com- 
pared with a Standard Map trajectory: the 25-KP trajectory 
spends considerable time trapped in a cell before escaping 
onto the next; the Standard Map trajectory looks like a sim- 
ple random walk. 

result is not fully explained. Here, we find D^o ^ K'^e in 
this regime. 

The structure of the paper is as follows: 
(a) In Section II we review the classical dynamics of the 
Standard Map and the 25-KP, we explain the cellular 
structure of the 2(5-KP, we show it is natural to work in 
rescaled momenta Pe — pe to remove effects dependent 
on the cell size 27r/e upon which there is an effective kick 
strength if ^ = Ke and we obtain an approximate map for 
the trapping regions and present a heuristic justification 
for the Doo ~ K^e behavior which is seen in the Ke < 1 
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FIG. 2: Classical diffusion of the Standard Map vs the 2S-KP. 
Here A" = 7 in both cases. For the Standard Map (broken 
lines) for K ^ 1, the same diffusion rate ~ /2 charac- 
terizes all timescales and all starting conditions. For the 25- 
KP (e = 0.05), the diffusion rate D{po,t) at first depends on 
the initial momentum po and time. The inset magnifies the 
< 50 behavior. For atoms prepared near the cell bound- 
aries, ie Po = ptrap ~ (2m + l)7r/e, where m — 0, 1,2..., we 
find D{po — Ptrap, t) — for the first 20 or so kicks, but then 
diffusion speeds up rapidly. For atoms prepared at the center 
of the cells, po — 2mn, there is initially rapid diffusion un- 
til the ensemble reaches the cell boundaries. We see that as 
t oo, (in this case, for N > 200) the diffusion rate reaches 
an asymptotic value D{po,t) = Doo independent of po . But 
Doo < Do: the diffusion rate is always slower than the uncor- 
related Standard Map. The quantum behavior is illustrated 
by the dotted lines. For effective h = 1/4, the quantum diffu- 
sion follows the classical behavior for only a few kicks before it 
is arrested by Dynamical Localization. For smaller h ~ 1/30 
it follows it for considerably longer. 



regime (where fractional power-law behavior was found 
for the quantum localization lengths L ^ h^^ "^^ [TlLll^ ). 

(b) In Section III we explain how to derive the form of the 
long-range correlations which modify the diffusion rate 
of the 2(5-KP. While these results were used to analyze 
the experiments in |lOj |. due to space constraints a full 
derivation could not be presented there. The emphasis 
is on timescales comparable to the quantum break-time 
t*: since the diffusive behavior is 'frozen-in' at this point, 
even transient (but long-lived) classical correlations be- 
come important for the quantum dynamics and hence are 
essential to understand the experimental data. This sec- 
tion is the most technical however, and readers without 
a detailed interest in classical correlations can obtain the 
key points from Figs. 5-7 and their captions. 

(c) In Section IV we investigate for the first time the 
asymptotic linear diffusion regime {t —^ oo). We attempt 
to derive analytical forms for Doo- For large Ke, one 
of the main results of the present work is the formula 

Doa — K^[l — Y^"/^^^] which gives extremely accurate 
results compared with numerics (see Fig|5}. It is shown 
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from this that the diffusion rate of the 25-KP can ap- 
proach, but not exceed the uncorrelated rate (ie the value 
expected if the particles executed a random walk). We 
attribute this to the absence of accelerator modes 0| 
which have been observed experimentally for the usual 
QKP 15] but are not found in the 2S-KP. For Ke < 1, 
although we show that including higher-order families of 
long-ranged correlations makes the analytical diffusion 
tend to the numerical result we are unable to quantita- 
tively reproduce the D^^ ~ K^e behavior, (d) Finally, 
in Section V, we conclude. 



II. CLASSICAL DYNAMICS 
A. Standard Map 

The classical Standard Map is obtained by integrat- 
ing Hamilton's equations for motion in the J-kick poten- 
tial V{x,t) — —KcosxJ2n^(^ ~ ^T). One obtains two 
equations which may be solved iteratively to evolve the 
system through each period T: 

PN+i ^ PN + K sinxN; XN+i = XN + Pn+iT 

(2.1) 

A 5-kick is followed by a period of free evolution with 
constant momentum. With increasing kick-strength K 
the system makes a transition from integrable, regular 
dynamics to eventual full chaos. For K > 1, chaotic 
momentum diffusion is unbounded and all chaotic phase- 
space regions are connected. If we neglect all correlations 
between impulses, ie assuming < sm.XNSm.XN> >~ for 
all kicks, the momentum of a trajectory in effect rep- 
resents a 'random walk'. The corresponding energy of 
an ensemble of particles grows linearly with time, since 

< p2 >^ j^^jy ^ K^/2N. NB: from Sec.III.B onwards, 
time is measured in kick-pairs t = N/2; in that case 
all diffusion rates are doubled; eg < p^ >^ Dot = K'^t; 
hence the uncorrelated diffusion rate in those units would 
be given by Dq — K^. 

The overall diffusion in the chaotic regime is in general 
however not uncorrelated (provided K is not too large). 
In the effect of short-range correlations between 

kicks was investigated theoretically. A more accurate 
form for the diffusion rate D ~ K^{^- J2{K)- jf{K)...) 
was obtained, where Jm{x) is a regular Bessel function 
of the first kind of order m and argument x. The sec- 
ond term is a 2-kick correction resulting from correlations 

< sinx^f sina;Af+2 >; the third term is a 3-kick correction 
resulting from < sinx^f sin a; +3 >. The effects of these 
corrections on the energy absorbed by atoms in pulsed 
optical lattices have been experimentally observed [l5l |. 
Note that for the Standard Map, the correlations repre- 
sent a simple change in the magnitude of D\ the energy 
increase is still linear in time. In jl6l | it was further shown 
that in an asymmetric potential, the 2-kick correlations 
yield a local correction to the diffusion, ie D depends on 



both time and the relative initial momentum, poj between 
the atoms and the standing wave of light. This produces 
a type of chaotic Hamiltonian ratchet. 

B. The 2(5-KP 

The classical map for the 25-KP is a straightforward 
extension of the Standard Map: 

PN+i =Pn + KsinxN] Pn+2 ^Pn+i + KsmxN+i 

XN+1 = Xn + PN+ie; Xn+2 = Xn+1 + PN+2T 

(2.2) 

where e is a very short time interval between two kicks in 
a pair and r is a much longer time interval between the 
pairs. It is easily seen from the map that atoms for which 
p^e — (2m-|-l)7r and m = 0, 1, 2, ... experience an impulse 
Ksvuxq followed by an impulse ~ Ksm[xQ + tt) which 
in effect cancels the first. The regime po — (2m. -I- l)'''/^ 
corresponds to the 'momentum trapping' regions. Con- 
versely in the case p^e = 2m7r, a series of near-identical 
kicks produces initially rapid energy growth. 

Some of the characteristics of the diffusion can be an- 
alyzed by the properties of the classical map in the trap- 
ping regions. Starting from the map 12. 2|) with = we 
re-scale all variables p'^ = pe and K^^ — Ke, — r/e ^ 1 
to obtain 

Pi=Po + sin zo; P2 = Pi + ^£ sin xi 

Xl^Xo+p\] X2^Xl+plTe (2.3) 

We take Pq — Pfi + Sp"^, where is the trapping mo- 
mentum, ie = (2to -t- l)7r and we choose m = for 
the first trapping region. Inserting the above into 12.31 we 
have 

P2 — Pa ^ sinxo -I- sin(a;o + + 5p'^ + sinxp) 
= Po + sinxo -I- Ke[sm[xo + tt) cos{5p'^ + sinxo) 
-I- cos(xo + tt) sm{5p'^ + sinxo)] 

(2.4) 

Assuming small angle identities throughout, cos(/(e)) ~ 
1 and sin(/(e)) ~/(e), 

P2—PQ — Ke cos xq [Ke sin a;o + ^p'^] 
K^ 

= Pq — KeSp'^ COS Xq (2.5) 

So at the center of the trapping region {Sp'^ — 0) the 
double- kick map is equivalent to an effective sin 2x single 
kick Standard Map: 

pf ~pI- KRsm2xo (2.6) 

where Kj^ — K^/2. Further away from the exact trap- 
ping momentum where Sp"^ » Ke/2 we have a cosinu- 
soidal map: 

pf - Pq- cos Xf) (2.7) 
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where Kgp = K^Sp^. 

From these arguments we clearly see that the natural 
stochasticity parameter of the 25-KP in re-scaled mo- 
menta p'^ = pe is iff = Ke. It is also important to 
estimate the range of < \5p'^\ < Sp^^^: the trap- 
ping regions have a small but finite width determined by 
whether there is significant cancelation between consecu- 
tive kicks. Requiring K^^ sin xq + sin(j:o -l- tt -|- 5p'^) ~ 
we estimate 5p^^,-^ ~ 7r/6, in other words over about a 
sixth of the width of each momentum cell, classical trajec- 
tories experience significant trapping in the momentum 
diffusion. 

In Figures 13 and 01 we can see the change from a sin 2a; 
map at p'^ ~ p|j where we observe two sets of stable is- 
lands within the range < x < 27r, to a cosinusoidal 
map further out where the position of islands is shifted 
by a phase of 7r/2. Fig|31 compares the detailed structure 
at p*^ ~ p'j^ with a standard map phase space for which 
KsM = Fig2] compares an enlarged section of a 

'side resonance' further away from p^, with a standard 
map with Ksm — KeSp. In both cases the correspon- 
dence is clearly visible. 
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FIG. 3: Correspondence between the local phase space at 
the resonant momentum pa of a trapping region {K — 4.8, 
e — 0.02, r = 1.98) in the double-kicked system (top) and 
a standard map (bottom) with a kick P2 ~ Po — -ffi? sin 2a;o , 
where Kr — = 0.24 is the kick strength of the Standard 
Map and the period is T = 2. Note that unsealed momenta 
are used here. 

We note that in the Standard Map there is a transition 
to unbounded diffusion for a critical kick strength K = 
Kcrit — 1 where the last invariant curve breaks, leaving 
a fractal cantorus structure. This curve corresponds to 
momenta p ~ 27rmM, where E is the golden ratio. In 
it was shown that for Kcrit < K < 4.5 the diffusion rate 
is D ~ 0.3{K — Kcrit)^- The value ~ 3 was obtained in 
[l^ from the scaling properties of the classical map near 
the golden-ratio cantorus. 

From the above, one can hazard a rough order-of- 




FIG. 4: Correspondence between the outer phase space near 
the resonant momentum pr of a trapping region {K = 4.0, 
e — 0.02, r = 1.98) in a double-kicked system (top) and a 
7r/2 phase-shifted cosinusoidal standard map (bottom) p2 — 
po — Ksp cos a;o , where Ksp = KeSp — . 75 is the kick strength 
of the Standard Map and the period is T = 2. Note that 
unsealed momenta are used here. 



magnitude justification for a cubic dependence of the 
asymptotic diffusion rate Doo — K^e for O.l < < 1 
found numerically (see Sec. IV). Note that in scaled mo- 
menta. Doc ^ K^- We note the correspondence between 
the Standard Map and the 25-KP with K ^ K^; for 
small K^, 5p^ > 0, the form of the map in 12.71 is domi- 
nant, since trajectories spend most of their time in the 
trapping regions. In this case, the scaled critical kick 
strength K^"* ~ e. Taking a representative value of 

< Sp" >~ (5p^„j./2 = 7r/12 ~ 1/4, the effective kick 
strength in 12. 71 is ~ K^/A. We thus have a diffusion rate 
in the trapping regions in terms of re-scaled momenta 

< >~ 0.3(-^ — K^'^^*^)^t. Then, in unsealed momenta, 
we would have D ^ K^e provided K is reasonably large 
compared with e. While this is by no means a rigorous 
argument it may provide some indication of the mecha- 
nism underpinning the dependence. We note the 0.75 
exponent in the L rx h~^ '^^ scaling is also an exponent 
associated with scaling properties of the region around 
the Golden ratio cantori. We propose tentatively this 
supports the proposal in 0, that the properties of 
the 2(5-KP are intimately connected with the properties 
of the Golden ratio cantorus. 



III. MOMENTUM DIFFUSION 

A. Standard Map 

The classical diffusion corrections for the Standard 
Map were first obtained by Rechester and White 
and we follow their notation closely. From the map 
(|2.1() . the momenta and positions of a trajectory evolve 
by a sequence of impulses: pn — Po + Sn-i and xn — 
xjv_i +po + Sn-1, where Si = X]j=o Ksinxj, the initial 
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momentum of an atom is po and the period T is taken turn G{xo,po, t = 0), at a later time (measured in number 
to be unity. If we consider an ensemble of particles with of kicks t = N) the distribution is given by 
an initial probability distribution in position and momen- 



+00 +00 „27r ^2ir 

G{xt,pt,t) = X! •■■ X! / dxodpo G{xo,Po,0) dxt ■■■ j dxi S{pt - po - St-i) 

X 5{xt - Xt-1 ~Pq - St-1 + 2nnt) ... 5{xi - xq - pa - Sq + 27rni) 



27r 



(3.1) 



The sums over ni ... appear because of the periodicity By taking the initial distribution as G{xo,po,0) 
of phase space in ... Xt. The momentum diffusion rate ^S{p — po) (ie a uniform spatial distribution with all 



D is given by 



D = hm - <{pt- pof >t 



t—^oo t 



dxt 



+ 00 



dptG{xt,pt,t)ipt-pof (3.2) 



particles at initial non-zero momentum pq) and using the 
Poisson sum formula giving the Fourier transform of a S- 
spectrum, J^n '^(2/ + 2nn) — ^ J2m e""^, we can rewrite 
El as 



1 



00 



D = lim / ... / 

t—>-oo t ^ ^ 

mf — ~oo mi— — 00 I— 



n 



dxj 
2^ 



{St-if exp{^imj(xj - x^-i - po - 'S'^-i)} 
i=i 



(3.3) 



To lowest order one can set all rrij coefRcients to zero, 
thus eliminating all exponentials. By using the previous 
form of S and integrating over the sine products it is 
easily seen that the random walk Dq = K"^ /2 is recov- 
ered. Higher-order corrections to the diffusion rate are 
obtained by setting certain rrij coefficients to a non-zero 
value; for the most dominant corrections \mj \ ~ 1,2. The 
integrals are solved using the relation cxp{±ii4rsinx} = 
Sri^oo Jn[K) exp{±ma::} and for corrections to be non- 
zero, all arguments of exponentials must vanish for 27r- 
periodic integration. This requires pairing exponentials 
with the relevant Bessel functions and sine products in- 
cluded in S. 

The main corrections to the Standard Map are the 2- 
kick and 3-kick correlations found in |l3l | and account 
in larg e measure for the experimental oscillations seen 
in [13 • The 2-kick correlation is obtained from set- 
ting rrij = ±1 and rrij-i = =f1 and the lowest or- 
der 3-kick correlation from ruj = ±1 and = Tl- 
From the sinxj suiXj-2 and sinxj sina;j_3 terms we ob- 
tain C2 = -K^J2{K) and C3 = -K^JfiK) respectively, 
as previously. There is also a higher-order 3-kick cor- 
relation, C3 = +K'^J^{K), found in derived from 
rUj = ±1, rrij^i — =p2, mj_2 = ±1, and a 4-kick cor- 
relation, C4 — +K^J2{K), cited in [23| derived from 
rrij = ±1, rrij^i = rnj-2 = tIj 'Tij_3 = ±1. This leads 



to a total correction to the diffusion of 

D = _ J2(X) _ Jl(K) + Jl{K) + 4{K)] (3.4) 

These terms represent the correlations between two given 
impulses only, two, three or four kicks apart. In the 
next section we shall see that in the double-kicked system 
there are entire families of terms representing correlations 
between a given impulse sina;^ and every other impulse. 
Such a global correlation family does in fact exist for the 
Standard Map, as the above method of derivation can 
be extended to any order of k-kick correlations between 
any two given kicks {smxj sinxj-k). For the lowest 
order we have rrij = ±1 and rrij-k+i = Tl and hence 
Cfc = -K^Jl{K)J^-^{K) (fc > 3). Clearly corrections 
become smaller with increasing k. The 2-kick correla- 
tion is a special case of this global collection of terms. 
Note the linear time dependence of all individual k-kick 
correlations (ie < p^ >cx Ckt). The total correction to 
the momentum diffusion rate due to all the above k-kick 
terms is Er=3 <^fe = -K'^iAiK)) / {I - Jo{K)). Similar 
higher-order global corrections can be found in the Stan- 
dard Map of the general form ±K'^ {]\ {K)) / - 
Jo(-ftr))") for some m,n,p. However in practice, higher- 
order corrections beyond the basic terms in can gen- 
erally be neglected for all K > 5 since they do not alter 
the diffusion rates significantly. For smaller K there can 
be significant differences; however phase space becomes 
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increasingly regular as K decreases and for < 2 a diffu- 
sive approacli is not justified. In the double-kick system 
the higher order (long-range) correlations originate only 
in the thin trapping layers so a diffusive analysis is effec- 
tive even in a regime where long-range correlations are 
important. 

It should be noted that none of these corrections de- 
pend on the momenta of the atoms in the ensemble; they 
simply alter the magnitude of the overall linear rate of 
energy absorption for a given K. In any (5-kicked system 
where the kicking periods are all equal and there are no 
other asymmetries, effects associated with momentum- 
dependent diffusion corrections are negligible. This is 
due to such corrections including oscillations cosTpo on 
a comparable scale to the natural width of the initial 
momentum distribution Apo ~ 2tt of the atomic ensem- 
ble. Hence these corrections average to zero. The 1-kick 
correlation sinxj sinxj-i between consecutive kicks for 
which a single rrij = ±1 only, is such a correction and 
is hence absent in the Standard Map, but will be seen 
to contribute significantly to the double-kicked system, 
where a much shorter kicking period is introduced. 



B. The 25-KP 



For the 25-KP the notation of the Standard Map 
diffusion equation can be changed slightly to include 
two closely-spaced kicks for each time step denoted by 
(1) and (2); thus the evolution of momentum is now 
in terms of pairs of kicks. Throughout this section 
we work in unsealed momenta < >= Dt (in re- 
scaled momenta < pi >— De^t). S takes on the form 

^i'^^ — J2ln=iJ2l=i^^^^^"^ E3l ^'^d indicate ex- 
plicitly the two time intervals r and e, defined as previ- 
ously. Again, Dq is obtained by setting all rrij coefficients 
to zero, but it should be noted that as there are now 
2t variables, the form of Dq changes to and hence 
< {p — pq)'^ >— Dot = K^t (or in scaled momenta K^t). 
Obviously this does not change the underlying physics; 
the new formula is only due to a redefinition of time in 
terms of number of pairs of kicks, t = N/2, for even N. 
Physical time is t{T + e). 



D = 



oo oo oc 



°° * r^TT , (2) ,27r , (1) 



2tt 



2tt 



X exp{J2[tmf\xf^ - xf^ - e{po + ')) + tmf\xf^ - xf\ - r{po + S^%))]} 



(3.5) 



C. Momentum-dependent Diffusion of the 25-KP 

The introduction of the short timescale e results in a 
significant 1-kick contribution to the diffusion, unlike in 

(2) 

any previously studied i5-kicked system. For rrij = ±1 
the correlation involves cospoe and Bessel functions of 
argument Ke — K^, while for mj"'^'' = ±1 the correla- 
tion involves cospo''" ^^nd Bessel functions of argument 
Kt. The latter case gives negligible contributions as r 
is the large time interval between pairs of kicks, result- 
ing in fast oscillations with pq. We can effectively set all 
m^-^^ = for all momentum-dependent correlations - a 

(2) 

valid approximation provided r >> e. The rrij coeffi- 
cients will be simply referred to as rrij for the remainder 
of the paper. 

For the case of cospoc, the sina;^^'' sinx^^^ term results 
in a kick-to-kick correction 

Cit = cospoe(Jo(i^e) - J2{K,)) 
t 

X J24^-\K,) (3.6) 



where t — N/2 {N even) ie time is measured in pairs of 
kicks as for all correlations that follow. Kicks at longer 
times have weaker correlations than those at short times, 
however since Jq{K^) ^ 1 as e ^ 0, the time-dependent 
summation decays slowly with time and the lifetime of 
the overall kick-to-kick correlation may far exceed the 
active running time of an experiment. Note however that 
the correlation does decay to zero eventually (as for all 
other momentum-dependent correlations), ie D{po,t 
oo) = 0. Hence we do not take the infinite time limit in 
13.51 instead we calculate the quantity Ct. The geometric 
sum in l3.6l saturates to 1/(1 — Jq) after a time ^ 10/{K^)'^ 
and hence 

hm C,t ^ X^cospo /°f^\~/;f^^ (3.7) 

For short times Cit grows linearly with time, since for 
small we have Jo{Ke) ~ 1 >> J2{K^). We can ap- 
proximate the correlation to K^tcospoe and hence the 
average energy of the double-kicked particle also grows 
linearly in this regime. 

Figgis the same as in [ifl] showing a numerical simu- 
lation of the energy absorption of an ensemble of 100,000 



7 



classical particles at K = 7, e = 0.05, t = T — e = 1.95, 
as a function of their initial momenta po various times 
(measured in pairs of kicks). The numerics are super- 
posed with combinations of the time-dependent corre- 
lations Cjt presented in this paper. In Fig[SJi, the ba- 
sic cosine behavior of Ci is clearly visible: atoms with 
initial momenta po — 2m7r/e = 0,125.66,... experience 
the largest energy absorption, while those prepared at 
Po ~ (2m -I- l)7r/e = 62.83,188.50,... are 'trapped' near 
this initial momentum and experience almost no energy 
absorption. When one looks at Figures and|Sli, some- 
thing unexpected occurs. The maxima of Fig[Sli slowly 
turn into near minima, while energy absorption for atoms 
near the 'trapping' regions increases continuously. A 
complete reversal of the initial situation eventually oc- 
curs at longer times, < {p — p^Y >oc — cospQe. 



The reasons for this lie in considering a whole fam- 
ily of global correlations similar to the case of the Stan- 
dard Map. The lowest order global correction originates 
from ?7ii = ±1, just as for C^, but this time we look 

(2) . (r) 

at sinx^ sina;^ terms, where r = 1,2 and i < j but 
otherwise arbitrary. In this way we include all k-kick 
correlations for fc > 2, ie correlations between the second 



impulse in a given pair, sinx^ 
pairs. 



(2) 



, and impulses in all other 
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Cat = ^ Cfet = -2K^ cospoeJ^K, 
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FIG. 5: Agreement between analytical difFusion corrections 
and numerical double-kick simulation of 100,000 classical par- 
ticles a.t K = 7, e = 0.05, r = 1.95, for various times mea- 
sured in pairs of kicks, t — N/2 where is the (even) number 
of individual kicks. Physical time is t{T + e) where e is the 
short time interval between kicks in a pair and r is the long 
time interval between pairs. Energy absorption is plotted as 
a function of initial momentum of the particles. The correc- 
tions included in the analytical curve in (a) are Cit, Ccit 
and the lowest order C^lt term labeled Ila in Table 1. In (b) 
all correlations given in Table 1 are included. Agreement is 
excellent at short times, but higher-order terms are needed at 
later times. One observes a sign reversal of the cosine enve- 
lope and gradual disappearance of the inverted Poisson peaks 
at initial minima. Particles initially prepared in momentum- 
trapping regions eventually absorb the most energy. At very 
long times the lowest order global correlation Cgi dominates 
almost completely. 



xY,i2j~2)4^-'{K,) (3.8) 



The corrected energy is now < (p — po)^ >— K^t + Cit + 
Ccit. Every individual correlation between two kicks in 
Cgi is small compared to the nearest neighbor correla- 
tion Ci, as Jf{Ke) « Jo{Kf). Importantly however, 
at a given time the correlations do not get weaker with 
increasing k but are equal in size for any non-nearest 
neighbor correlation. All correlations however do be- 
come weaker with increasing time and eventually sat- 
urate. Summing over t reference impulses and 2j — 2 
paired correlations for each, results in a total correction 
to the diffusion which after a given time becomes domi- 
nant relative to Ci . Notice the difference in sign between 
the 1-kick correlation and all others. The summation in 
13.81 can be seen to be the derivative of that in 13.61 and 
hence 



Cat = -AK'^cospoeJ^iK,) 



(3.9) 



For small it can be shown by assuming Jo{Ke) ~ 
1 — and expanding the above to second order binomial 
that Ccit initially increases quadratically at short times: 
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l~K,~t{l~ + {t- 1)(1 - is:,)2*+i ~ 1 - if, - <[1 - (2t - 1)K, + (2f -3t+ 1)K^] 

+ {t - 1)[1 - {2t + 1)K, + (2t2 + t)K^] = 2K^{t^ - t) (3.10) 



For the parameters in Fig[Sl Ccit saturates to a value 
approximately twice that of Cit as i — > oo. 

hm Cait = ^AK^cosp,e^^^^^f^ (3.11) 

The behavior of the 1-kick and all global correlations 
with time is illustrated in FigEl (including higher-order 
Poisson correlations presented later). The absolute 
value (|C|t) of the maximum energy absorption (where 
cospo£ — 1) with increasing number of kick pairs is 
shown. Note that all momentum-dependent correlations 
have a non- linear time dependence and saturate after suf- 
ficient time. This means that the momentum-dependent 
diffusion rate approaches zero at long times, where all 
particles absorb the same amount of energy irrespective 
of pq. At short times the effect of Cgi is negligible com- 
pared to Ci, but at later times the global correlations 
dominate the diffusive process, explaining the sign re- 
versal in Fig[3i. The importance of the global correla- 
tions persists for values of for which the correspond- 
ing phase space is completely chaotic, a phenomenon not 
previously observed. 

FigEb corresponds to a regime where Ci and Cqi are 
of similar importance (near the crossing point in FiglSJ. 
Here and in Figl^J: another feature mentioned earlier be- 
comes evident. The initial troughs of the cosine in FiglSja, 
turn into narrow downward peaks, superimposed onto 

I 



where /(j) depends on the particular term: some exam- 
ples of Poisson terms are given in the appendix. Each 
additional mu = ±1 coefficient adds two multiplicative 
Ji factors of possibly different arguments to both Poisson 
families, since Bessel function arguments can now be any 
integer multiple of i^e, depending on the exact combina- 
tion of rrik coefficients. This continuously decreases the 
significance of the correlations, since Ji{x) — > as a; ^ 0. 
Hence the dominant corrections are those of lowest order 
in Ji (and with few non-zero mj. coefficients). Note that 



I 

the cosine envelope. The origins of these peaks are global 
correlation terms of higher cosine orders cosnpoe, which 
yield near (5-functions at the relevant momenta through 
X;„(-l)"cosnpoe = Em ^(Poe - (2to + l)7r). Such Pois- 
son terms arise when more than one rrij coefficient in 13.51 
is set to ±1, and Ej "^j = the sum of the coefficients 
defines the cosine order n. Terms with \mj\ > 1 include 
Bessel functions of increasingly higher orders; since for 
e ^ 0, Jn{Kg) more rapidly for increasing n > 0, 
such terms can be neglected here. 

Derivations of higher-order global correlations are te- 
dious, yet reasonably straightforward. There are two 
distinct families, one following the pattern of Ci and 
the other following Cqi, each causing different diffusive 
behavior. In both cases the terms of interest involve 

(2) . (r) 

sinx^- sin a; J. where r = 1,2 and k < j. In the for- 

(2) 

mer case m^, — ±1 and correlations involve a Jq — J2 

(21 

factor as for Ci. In the latter case mj, = and cor- 
relations involve a 2Ji factor as for Cgi- The behavior 
depends on whether the correlation is with a kick in a pair 
associated with a corresponding zero or non-zero co- 
efficient. We choose to denote the global Poisson terms 
as CqI and C^l depending on the family they belong 
to, where n is the cosine order. Cgi is the lowest order 
global correlation family belonging to Cq^. All Poisson 
terms have the following general forms for some m 



(3.12) 



(3.13) 

I 

for Poisson terms other than Cgi various k-kick correla- 
tions are of similar but not necessarily equal strength at 
a given time t. 

For n — I there are higher-order Poisson terms (in Ji) 
contributing only to the cosine envelope in Fig|31 These 
are shown in Fig|H|i in comparison to the lowest order 
corrections Cit and Cgi^- Again absolute values of max- 
imum energy absorption are shown. As 0{Ji) increases, 

terms become less significant. For all terms, Cqi cor- 
relations eventually become more dominant than CqI 
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FIG. 6: Behavior of global correlations Ccnt with number of 
pairs of kicks. Absolute values \C\t of maximum energy ab- 
sorption are shown (cosnpoe = 1). FiglHK shows the compari- 
son between individual correlations at order n — 1, including 
nearest neighbor Ci, lowest order global Cci and higher order 
Poisson correlations. FiglSJj shows the comparison between 
different cosine orders n. Trapping terms C^^ are shown with 
stars, absorption enhancing terms Cq^ with circles. The la- 
bels correspond directly to those in Table 1 showing explicit 
values at t = 15. Terms are denoted by the cosine order n 
resulting from the total sum of non-zero rrij coefficients, and 
the order in Ji resulting from the number of non-zero coeffi- 
cients. Curves ininji are made up of the various correlations 
in Table 1 for a given order n. Note the linear and quadratic 
rise respectively for Ci and Cgi at early times. Absorption 
enhancing terms always overtake their momentum trapping 
partners at equal order after sufficient time, but the crossing 
point shifts to later times as both the order n and 0(Ji) in- 
creases. Overall importance of terms also decreases for higher 
orders. 

correlations, but the crossing point shifts to later times 
with increasing 0(Ji). The difference in saturation ener- 
gies between the two types of correlations also becomes 
smaller. 

FigEb shows the comparison between the combined 
n = 1 correlation families and higher cosine orders. The 
behavior with increasing n is as before; correlations 
always become more dominant than their partners of 
equal order, after sufficient time. Cq^ terms increase 

as ~ t", while Cq^ terms increase as ^ 

Table 1 shows the maximum energy absorption values 
of the dominant Poisson terms at t = 15 (30 kicks) cor- 
responding to Figlsjs. The most significant diffusion cor- 
rections are those for small n and 0(Ji) and importance 
of terms rapidly decreases with higher orders. Terms of 
0{Jl'^) are about 50 times smaller than the leading cor- 
rections. 

The momentum diffusion corrections derived here now 
enable us to explain the behavior seen in Fig|Sl including 
the inverted peaks. The agreement between the numeri- 
cal data and the analytical diffusion calculations is gen- 



9 



Term 


rtij pattern 


0(Ji) 


Value (t = 15) 


la (Ci,Cgi) 


±1 


0,2 


-f472,-355 


lb (4\',4i') 


±l,±l,=Fl 


4,6 


+100,-33 




±l,±l,±l,Tl,Tl 


8,10 


-1-14,-4 


Ila {Cq2 ,Cq2 ) 


±1,±1 


2,4 


-227,-1-113 


lib (4i,4i) 


±l,±l,±l,=Fl 


6,8 


-59,+24 


Ilia (Cq3 ,C(J3 ) 


±1,±1,±1 


4,6 


+82,-39 


iiib (4i,4i) 


±l,±l,±l,±l,Tl 


8,10 


+29,-12 


IVa (4i,4i) 


±1,±1,±1,±1 


6,8 


-29,+14 



TABLE I: Diffusion correlations shown in FiglH] 



erally good. At very short times (FiglS^t) essentially only 

Ci, Cqi and Cq\ (Ila) contribute to the diffusion and 
agreement is excellent. In Fig|3lD all of the correlations 
in Table 1 have been included in the analytical curve 
and good agreement is obtained: the sign of the cosine 
envelope starts to change and the inverted peaks slowly 
vanish at later times. From ITT^ and IXTSl one notes that 
at the same order n, the two types of correlations are 
of different sign and so oppose each other, clearly seen 
in the case of n = 1. For n > 1 \t is found that while 
Cq\ correlations contribute to increasing the size of the 
downward peaks in all cases and thus favor momentum 
trapping, C^^ correlations result in enhanced energy ab- 
sorption for initially trapped atoms. These latter corre- 
lations cause the release of atoms from trapping regions 
in momentum space. 

In FiglSJ; the analytical curve has been omitted as rea- 
sonable agreement cannot be achieved using only the cor- 
relations in Table 1. Higher-order terms are needed at 
intermediate times, however at very long times the dif- 
fusion is dominated by Cqi, the most important correla- 
tion, as can also be seen in FigE^L. After saturation, the 
higher-order n > 1 Poisson correlations, result in super- 
imposed upward peaks at the maxima of the cosine enve- 
lope giving the latter a pointed appearance. The newly 
found global correlation families lead to a situation where 
at long times, when diffusion no longer depends on ini- 
tial momentum, those atoms that started in a trapping 
region have actually gained more energy than those that 
started in an enhanced absorption region in momentum 
space. 

So far this paper has only treated the diffusion prob- 
lem classically, but a real experiment would obviously be 
carried out in the quantum regime. In ^3 ^^'^ 
it was shown that the effects presented here can read- 
ily be observed in a quantum experiment, and indeed a 
range of other interesting features were discovered. 

As mentioned in the introduction to this paper, energy 
absorption does not continue indefinitely in the quantum 
case. The energy saturates to a near-constant value after 
a characteristic quantum break time t* ~ , where 

?i is a scaled Planck constant. The smaller fi, the longer 
the system follows the classical predictions (see Fig|21). 

Fig[7| shows the energy absorption of a cloud of cold 
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atoms for K = 3.3 and fi, = 1, measured in the 25-KP 
experiment in Values of e vary from 0.045 to 0.160 
such that = 0.1485,0.3102,0.528. The resemblance 
to Fig[Slis evident and is due to these different values of 
Kf^ which control the relative importance of diffusion cor- 
relations. It should be noted that in experiments, e and 
r are usually dimensionless scaled quantities, such that 
e = A/T, where A is the physical time between kicks in 
a pair and T is the physical period of the double-kick sys- 
tem. The measurements are all taken well after the break 
time t* ~ Dq ~ 30, when the energies have saturated and 
no further evolution takes place. The experimental be- 
havior then depends on which diffusion corrections are 
dominant at the break time. As i^e decreases the time 
needed for the global diffusion corrections to dominate 
over Ci increases, hence in FigEK the system is arrested 
at a time where Ci is still hugely dominant. In Figures 
[7|d and[7|: however, is larger and the global correla- 
tions Ccn become more important when the break time 
is reached. In these cases the Poisson peaks and sign 
reversal for the cosine envelope are clearly seen. Fig- 
ures [7li-[7f show Ci and Cqi for the parameters in a-c, 
and the crossing point, indicating the point at which the 
global correlations become dominant, clearly shifts from 
to > t* to ta < t* at higher K^. Saturation values of 
correlations also decrease with increasing time between 
the double kicks, which results in the limit of the Stan- 
dard Map where e = t ~ 2 and there is no momentum 
dependent diffusion at all. Note that at very high values 
of K, and hence K^, correlations become less significant 
due to large Bessel arguments {Jn{x) ^ as x ^ oo). 

From the momentum distribution of the atomic en- 
semble one can also observe the trapping regions and 
the effects of the global diffusion corrections. After the 
quantum break time, the distribution of the usual i5-KP 
localizes exponentially in the momentum basis, causing 
the characteristic triangular shape on a logarithmic plot. 
For the 2S-KP the basic exponential shape has a 'stair- 
case' superposed onto it (see [TllIT^ '). 



IV. ASYMPTOTIC DIFFUSION: t oo 

Classically, the momentum-dependent diffusion in the 
2(5-KP does not continue indefinitely but is only a tran- 
sient effect at short and intermediate times. At very long 
times in the asymptotic regime, momentum-dependent 
correlations saturate and the diffusion rate is the same 
for all starting conditions po- In this regime diffusion 
is controlled by linear corrections independent of initial 
momentum, as for the Standard Map, which decrease in 
magnitude with increasing separation between kicks but 
remain constant in time. 

Such corrections are obtained from Eq. I|3.5I) in the case 
where the total sum of rrij coefficients is zero; hence n = 
and cos npoe =1. It is found that all these corrections 
are of opposite sign to the random walk Dq — and 
hence lower the overall rate of energy absorption for all 
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FIG. 7: Experimental realization of the 25-KP with a cloud 
of cold cesium atoms, pulsed periodically by pairs of laser 
kicks. The graphs show energy absorption of the ensemble 
for K = 3.3, e = 0.045, 0.094, 0.160 and h ^ I. Measurements 
are taken after the quantum break time, t* ~ 30, when the 
energies have saturated and momentum diffusion has been 
terminated. Varying Kf, means that the timescale to on which 
the global diffusion corrections become dominant varies in 
relation to t* . In FigUJi t* « to and the system is arrested 
in a regime where Ci dominates, corresponding to Fig|KJi.. 
Atoms prepared in the trapping regions remain trapped and 
absorb no energy. In FiglTJi t* ~ tc and the effects of the 
global corrections Con become visible: inverted Poisson peaks 
appear at trapping momenta and a partial sign reversal of the 
overall cosine envelope is seen. In Fig[71; finally, t* > to and 
the system follows classical energy trajectories long enough for 
the global corrections to become dominant. Trapped atoms 
start to absorb energy and diffusion in other parts of phase 
space approaches a (— cospoe) relationship. Figures d-f show 
the behavior of Ci (stars) and Cgi (circles) with time for the 
values of in a-c (r = 2 — e). We can consider to to be near 
the crossing point of the two correlations shown. This clearly 
shifts to earlier times as increases. 



values of K^. This is due to the absence of terms for 
which = 2 or m^j''^ ^ 0. Although for momentum- 

independent correlations, terms can also depend on the 
long timescale r between kick pairs, ie ti^-^'' = ±1, sim- 
ilar to the Standard Map case (r ~ T), such terms are 
far less significant than the dominant e-dependent correc- 
tions. Hence k-kick correlations for even k, such as the 
important 2-kick correction (oc J2), are effectively absent 
in the 2^-KP. The general form of correlations dependent 
on e only is 

t fU) 

m j oCfn—O m 

(4.1) 

Note the difference in sign in comparison to Dq and ab- 
sence of Bessel functions of order p > 1. Note further 
that although we denote these correlations by Ct, the 
actual diffusion rate D{t — > 00) is non-zero and can be 
found by taking the infinite time limit of 14.11 From the 



lowest order correlation for which 



±1, nik = Tl 
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(fc < j) we obtain a global family similar to the case 
of the Standard Map, including the usual lowest order 
3-kick correlation 

CGot{±l,Tl) = Jo2"(i^,) (4.2) 



=2 a=0 



and hence, for K^^l, 



D « K^[l 



(4.3) 



Note the similarity to the lowest order global family 
Ck presented for the Standard Map earlier. The mi- 
nor difference of a Jg factor stems from the fact that 
this family includes only the case of odd k. Higher or- 
der corrections for the 25-KP have the asymptotic form 
D = -K''Y{^{{Jl{mK,))/{l - J^{mK,))r for some 
TO, n. 

FiglSJi shows the effect of the two lowest order global 
correlation families on the energy absorption of the 
kicked particles for the same parameters as in Fig[51 
(ifj = 0.35). The change from the random walk 
is seen to be very significant. It should be noted that at 
early times Coot behaves non-linearly; however, this is 
only due to the inclusion of entire global families of corre- 
lations between all kicks. As time increases, more corre- 
lations are included in the formula, however each individ- 
ual k-kick correction still has a linear time-dependence, 
as for all correlations independent of momentum. At long 
times the energy increase becomes linear again, since cor- 
rections become weaker with increasing k. 

Fig|S|3 shows absolute values of the two corrections in- 
cluded in the analytical curve in in comparison to 
Do = ■ It is seen that the lowest-order term given in 
14.21 makes up the vast majority of the diffusion correc- 
tion. The second term (given in the appendix) is derived 
from a sequence of mj = ±l,±l,=pl,=pl and is hence 
0(Jf). The agreement between the analytical and nu- 
merical data in|HJi is reasonable and can be improved by 
including higher-order terms. 

Fig|5| shows a comparison of the ratio D/Dq between 
the Standard Map and the 2(5-KP. FiglSJ) shows the well- 
known diffusion behavior for the case of the Standard 
Map or (5-KP where Dq = and D is given by lO 

(this curve was first presented in |0 but without the 
C4 ~ K^J2iK) term). D/Dq vs K is shown and the be- 
havior is oscillatory around the usual random walk with 
some regimes of enhanced and some of hindered diffu- 
sion. At very large K the diffusion approaches the ran- 
dom walk value asymptotically. As was mentioned in 
Section II, higher order global corrections do not alter 
the shape of the curve in^p appreciably. For A' < 5 a 
linear relationship is found where D oc K'^, however as K 
is decreased the system becomes increasingly regular and 
the diffusion corrections no longer describe the behavior 
accurately. 

FigEli shows D/Do vs for the 25-KP, where Dq = 
(in unsealed momenta) and a series of analytical 
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FIG. 8: Effect of momentum-independent corrections to the 
overall rate of energy absorption by the 2(5-KP. Fig|H^ shows 
the significant decrease in energy absorption from the basic 
random walk. Do = . The analytical curve is composed 
of the corrections in Fig|H|D and shows reasonable agreement 
with the numerical data. Parameters are as in Fig|S] Absolute 
values of the correlations are plotted in|Hj3, so the analytical 
curve in|H^ is K^t - Caot{±l, =Fl) - CGot(±l, ±1, Tl, Tl)- 



curves are shown, including an increasing number of 
global correction families, in comparison to the numeri- 
cally obtained diffusion behavior. The analytical curves 
are denoted by the maximum order of Ji{mK^) included. 
The behavior for large is found to be different from 
the case of the Standard Map (seen in|^): for the 2d-KP 
the diffusion is always 'hindered' compared to the ran- 
dom walk; ie D < Dq for all K^, while for the Standard 
Map, there are maxima in the diffusion where D > Dq. 
For the latter, the maxima correspond to values of the 
stochasticity parameter K 2ttM where M is an integer. 
These coincide with the appearance of small transport- 
ing islands, (ie trajectories in these islands gain w 2tt 
in momentum for each kick), the so-called "accelerator 
modes" see eg the surrounding chaotic phase space 
generates Levy flights. This regime of anomalous diffu- 
sion was observed in the experiments in 15]. We have 
found no evidence whatsoever of any accelerator modes 
for K « 27rAf; we cannot exclude a small contribution 
from other higher-order accelerator modes, but the small 
maxima of the D/Dq in the 2(5-KP correspond to zeros 
of Ji {K^ ) and the diffusion is bounded from above by the 
random walk rate. Asymptotically, D — > Dq as — > 00 
as all correlations decay away completely. 

For small K^^ we again have a near-linear section where 
D oc K^e (or D (x in scaled momenta p"^) as men- 
tioned in Section II. The agreement between analytical 
and numerical results is excellent for K^^ > 2, where even 
the lowest order {0{Ji)) correlations are sufficient to re- 
produce the diffusion curve. As if — > agreement can 
be improved by including successively higher order cor- 
rections. 



12 




1.5 

1 

0.5 









(b) 


Standard Map 



10 



15 



20 



K 



FIG. 9: (a) The diffusion rate for the 25-KP (shown as a 
function of stochasticity parameter = Ke) never exceeds 
the 'random walk' rate Do; the uncorrelated rate, (for time 
counted in kick-pairs), is -Do ~ ■ The figure shows that 
for large Ke, there is excellent agreement between Eg 14.31 and 
numerics; in the regime of cold atom experiments, A'j < 1, 
dominated by fractal phase-space structures, agreement with 
the numerical result is improved by including higher order 
correlations (0( J")) but remains non-quantitative. Note that 
the numerical result D ~ is represented by a straight 

(dashed) line since D/Dq = Kt/2) 

(b) Illustrates the 'textbook' Standard Map diffusion rate 
D ~ K^/2[l- 2J2{K)...] obtained by Rechester and White 
fl^ . The rate oscillates about the 'random walk' value and 
the maxima occur where classical accelerator modes are im- 
portant. 



V. CONCLUSIONS 

We have analyzed the chaotic classical diffusion of the 
double-(5-KP. Although a straightforward extension of the 
Standard Map, this system exhibits a rich variety of new 
features. 

One motivation of the present work was to further un- 
derstand the behavior observed in experiments 10] and 
to p resent a more detailed derivation of the formulae in 
[lO'|. The experiment probed regimes where the classical 
diffusion is highly non-linear in time and depends on ini- 
tial conditions, but which were nevertheless sufficiently 
generic to be analysed in terms of corrections to a diffu- 
sive process. This is in contrast to the Standard Map, 
where long-range correlations are a feature of the near- 
integrable regime where most trajectories are stable so 



analysis based on diffusive processes is not useful. Al- 
though the non-linear regime of the 25-KP is, classically, 
transient, it is sufficiently long-lived 100 — 200 kicks) 
so that it is the only regime sampled by the quantum 
dynamics of the experiment in (IQj . 

Quantum studies of the 2<5-KP in the regime 0.1 < 
Kf^ < 1 were originally undertaken because they coincide 
roughly with the experimental values. However, the the- 
oretical studies then revealed novel quantum behaviour 
in the energy level-statistics and the fractional scalings 
L ^ /i"'^ ''^ of the quantum localization lengths J^J . 
It was suggested that the global quantum properties of the 
2(5-KP are, in this regime, intimately connected with the 
local scalings around the neighborhood of the Golden ra- 
tio cantori, as transport from cell to cell is limited by the 
former. The dependence of the diffusion rate found 
here lends further credence to this suggestion, as a cu- 
bic dependence is a feature of diffusion near Golden-ratio 
cantori |20j| . Unfortunately, as for the Standard map, the 
cubic behaviour cannot be derived analytically. At best, 
it was seen that adding successively higher orders of dif- 
fusive corrections makes the analysis tend towards the 
numerics. 

One of the main new results here is contained in Fig|51 
and Ea l4.3l in the regime ^ 1, where effects of can- 
tori and trapping are negligible. Here, the diffusion of 
the 2(5-KP follows a simple analytical expression but still 
differs strikingly from the Standard Map, in that it never 
exceeds the random walk rate (though it can equal it). 
The 2(5-KP has no two- kick correction (the J2{K) term 
which in the Standard Map roughly 'tracks' the accel- 
erator modes). We have found no numerical evidence 
for accelerator mode behavior. Hence, although there is 
much evidence for non-linear diffusion (of finite though 
prolonged duration) in the 2(S-KP, the major source of 
anomalous diffusion seen in the chaotic Standard Map is 
in fact absent; in this sense, 25-KP diffusion is closer to 
the random walk (as shown in Fig. 9) once trapping by 
cantori becomes unimportant. 

Finally, it is worth noting the potential applications 
of the 2d-KP, in particular as a velocity-selective atom 
filter in, for example, devices like an atom chip. Narrow 
trapping regions could be used to select atoms with 
Pq ~ 7r/e, while others would accelerate through the 
system nearly unperturbed. This could also be used to 
create very pure Bose-Einstein condensates, if located 
in a trapping region. A much stronger momentum- 
dependent effect is seen in the double-kick system than 
in the previously studied perturbed-period system |16| . 
Other applications in atomic manipulation may also be 
possible given further investigations. 

M Stocklin acknowledges a Graduate Research Scholar- 
ship; we acknowledge helpful discussions with Shmuel 
Fishman and Charles Creffield. 
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APPENDIX: GLOBAL CORRELATION FAMILIES 



Exact forms of two of the momentum-dependent global correlation families presented in this paper are given here 
as examples, as well as one further example of a higher order Cqo correlation term. 

(2) 

Note that the ' patterns given here and also shown in Table 1 correspond to the sum of correlations with 
such non-zero coefficients, including all possible permutations. The sum of the coefficients defines the cosine order 
n and permutations for which the partial sum, when coefficients are added from highest to lowest j, is zero, lead to 
zero-valued correlations. 

mj = ±1, ±1, =Fl (Term lb in Table 1) 



C^,t^-K cosp,eJ,{K,)J,{2K^)[Z j-^^ + 2 -j-^^ ] 



t i-3 



xEE-^o"(^^)^o^'"'""'(2X,) (A.l) 



±1, ±1 (Term Ila in Table 1) 



t j — '^ 

Jl^a^ Jo(Ae) Jo(2K,) 

— ±1, ±1, =f1, T1 (second higher-order correlation term in FiglHb and included in FiglHti) 



t J— 4 oi j-4— (Qi+a2) 

CGot{±h±hThTl)^-K'[JtiKe)JU2Ke)]J2J2 E E JS"HKe)JS"'{2K,)j',-^{K,) (A.3) 

j=4qi=0 02=0 "3=0 

Other Poisson and higher-order global correlation terms have similar forms and are derived in the same way as 
the above from 13.51 Note that the time behavior (C^^ terms increase as ~ t", while C^^ terms increase as ~ 
for small and t, can be shown for all terms by evaluating the sequence of geometric sums and expanding powers 
binomially to the appropriate order. Saturation values and asymptotic behavior are also straightforward to evaluate 
by taking the limit as i — > 00. 
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